The Carbon Flash 


Because of the strong electrostatic repulsion of carbon and oxygen, 
and because of the plasma cooling processes that take place in a 
degenerate carbon-oxygen core, it is extremely difficult to fuse car- 
bon. As long as the core remains under the Chandrasekhar limit, 
degeneracy pressure can support the star, and nuclear reactions will 
not occur. Stars with initial masses M < 7Mo, will probably stay 
under this limit, and thus their nuclear evolution will end with car- 
bon. 


If helium burning creates a core that is greater than the Chan- 
drasekhar mass, the core will collapse, and carbon fusing will occur. 
In objects with initial masses 8 $ M S 10Mo, this fusion will oc- 
cur under extremely degenerate conditions; the result is an explosive 
carbon flash. (The timescale for this burning will be milliseconds!) 
A deflagration front will probably develop, in which the shock wave 
itself does not ignite carbon, but instead, carbon is fused by the 
energy transported behind the shock front. (The alternative to this 
is a detonation front, where the shock itself ignites the fuel.) Either 
way, the star does not have time to adjust to the energy input, and 
the energy released (Q ~ 2.5 x 101" ergs gm7') is enough to disrupt 
the entire core. 
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Core evolution for high mass stars after helium-core burning. 


Carbon Burning 





In higher mass stars, the core temperatures are sufficient to fuse 
carbon under non-degenerate conditions. When this happens a 12C 
nucleus fuses with another !7C nucleus to form a compound excited 
state of 74Mg through one of many resonances that exist. The 74Mg 
then decays through one of many channels: 








Ce 2 a Mey Q = 13.93 MeV 
—» 23Na + !H Q = 2.238 MeV 

—» 20Ne + 4He Q = 4.616 MeV 

—> BMg +n Q = —2.60 MeV 

—, 160 + 4He + 4He Q = —0.11 MeV 

(26.1.1) 


with the electromagnetic and 3-body channels having the lowest 
probabilities, and the neutron reaction requiring the highest particle 
energy. Thus, to first order, '*C fuses to ?3Na and 7?Ne. However 
the protons and a particles quickly fuse as well, with many possible 
combinations. ‘The most likely are 


12C +H — BN +y 
BN =e Ee E 


13C + tHe — 1O +n (26.1.2) 


(This last reaction is important, since it creates a free neutron, 
which can be used to synthesize heavier elements.) In addition, the 
a particles may also interact with '7C, 160, 2°Ne, and ?°4Mg. The 
net result is that even though !*C fusion itself does not produce a 





large amount of energy, the additional reactions raise the amount 
of energy liberated to ~ 13 MeV per reaction. 


Although carbon and oxygen are roughly in equal abundance in the 
core, reactions between !*C and !°O are not important, due to the 
additional coulomb barrier. (By the time the temperature is high 
enough to fuse carbon with oxygen, all the carbon is gone.) Thus, 
the next series of reactions that take place are resonance reactions 
involving two oxygen nuclei: 











169 + 169 —, 32254 Q = 16.539 MeV 
—, 3IPh + 1H Q = 7.676 MeV 
—+ 31gj 4. Q = 1.459 MeV 
—, 289; 4 4He Q = 9.593 MeV 
—> ?4Me + tHe + “He Q = —0.393 MeV 


As in the case of carbon burning, the protons, neutrons, and a@ 
particles produced by this reaction can interact with a host of other 
nuclei, and thus a large reaction network must be analyzed. The 
primary product of oxygen burning is 7°Si, but a wide range of 
elements are produced. 


Photodisintegration 


At temperature of T ~ 10° K (which is typical of carbon and oxy- 
gen burning), nuclei can become disassociated by the thermal pho- 
ton bath. This is the exact analog of photo-ionization of electrons 
(except, of course, that the temperatures are very much higher). 


Photodisintegration rates can be calculated in the exact same man- 
ner as photoionization rates. Consider the reaction 


Kr a (26.2.1) 


In thermodynamic equilibrium, the photodisintegration rate, A,, 
must be related to the reaction rate Axy by 


NzaAy = Nx Ny Axy 


and the ratios of the species X, Y, and Z must be given by the 
Saha equation 


NxNy\ — GxGy (2mpkT\*? _ oper RT 
Nz rai 72 e ap 














where the G values are the nuclear partition functions, u is the 
reduced mass of the system, and @ is the binding energy. Thus we 
have 


CxGy (Inn I y A 
y = Gs ( a ) e OT yyy (26.2.3) 








GxG 11.605 
= 9.943 x 10°? “a, au exp (-= Axy 


where A is the reduced atomic weight, and Q is in MeV. 


During carbon-oxygen burning, 13N and ?°Ne are the principal vic- 
tims of photodisintegration. At temperature Tg > 0.75, the pho- 
todisintegration of tN is faster than its beta decay, and thus the 
concentration of 13C is sharply reduced. Similarly, at Tọ > 1.3, the 
photodisintegration of ?°Ne into oxygen and helium is greater than 
its creation rate from oxygen and helium. Thus 


Interestingly, the a particle produced by this reaction is now free to 
fuse with a different ?°Ne nucleus, giving 


20Ne + tHe —> *4Mg + 4.583 MeV 


Thus, ° Ne only exists for a short time in the star: it is made during 
carbon burning at T9 ~ 0.7 and destroyed during carbon burning 
at To ~ 1.3. 


Proton number 





6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 
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The nuclei which participate in the reaction network during silicon 
burning. The solid dots designate stable nuclei (i.e., observables); 
the open dots denote unstable nuclei. (See Truran, Cameron, & 
Gilbert 1966, Can. J. Phys., 44, 576.) 


Photodisintegration Rearrangement 


Since carbon and oxygen burning generate products such as 74Mg, 
28Si, and 3?Si, it is tempting to consider the next step where these 
elements get fused. However, the coulomb barrier for these species 
is so large, that photodisintegration will occur before their fusion 
temperature is reached. Instead, further elemental fusion occurs via 
protons, neutrons, and a particles that have been photodisintegra- 
tion from nuclei with small binding energies. 





Now consider the reaction specified in (26.2.1). If the photon bath is 
high, this reaction will attempt to come to equilibrium, and establish 
the abundance ratios given by the nuclear Saha equation (26.2.2). 
However, once removed from nucleus Z, particle X may also fuse 
with a different nucleus, P, i.e., 


PX —Q 


and nucleus Q may have a binding energy high enough to resist 
photodisintegration. Thus the abundance of element Q may be 
increased, and that of element Z may be decreased by photodisinte- 
gration rearrangement. This may, or may not create nuclear energy 
(since many of the reactions involve energy losses from neutrinos), 
but it will certainly relocate nucleons into those nuclei that are the 
most tightly bound. 


The results of photodisintegration rearrangement depend on the 
exact temperature and timescale of the burning process. Many 
of the species formed during the process will be radioactive, and 
their positron-decays will reduce the total proton-to-neutron ratio. 
(Again, the importance of these decays will depend on the timescale 
of the burning through their temperature dependence. ) 








The computation of nucleosynthesis during photodisintegration re- 
quires following a series of long differential equations containing all 
the reactions that can create or nuclei, t.e., 
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(Note that the above equation does not even include reactions such 
as X +p — Y +n, X +n — Y +a, etc.) Much of the work in 
setting up the network of equations is in deciding which reactions 
are important, and which can be neglected.) 


At oxygen burning temperatures, nuclei such as 32S, ?'Ph, and 2°Si 
are all destroyed in favor of ?SSi, which is the most tightly bound 
nucleus in the intermediate mass range. Once the temperature be- 
comes high enough to photodisintegrate 7°Si (To Z 3), a host of 
(a,7), (p, y), (n,y) reactions (and their inverses) begin simulta- 
neously. This is the Silicon Burning phase. Some of the fastest 








reactions run to equilibrium, but most reactions will not reach this 
stage. Thus, the results of silicon burning depend on how far to- 
wards completion the reactions run. During this time, there is also 
a slow leakage of nuclei from the intermediate-mass region to the 
iron group (which are the most tightly bound nuclei). 





As silicon burning nears completion, the nuclei come closer and 
closer to establishing nuclear statistical equilibrium. In this case, 
the nuclear Saha equation defines the abundance of each species. 
This has a very simple form. Consider that, from the Saha equation, 
the abundance ratio of species A—1, Z to A, Z is 


N(A-1,Z)Nn — GnG(A-1,Z) (2npkT\*°” Qı 
NAZ CAZ) = ) = (ar) 
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where G is the nuclear partition function (2 for the neutron), Q1 is 
the binding energy of A—1, Z, and © = (27m,kT/h?)3/2. Similarly, 


we can write an expression for the ratio of species A—2,Z—1 in 
terms of A—1, Z 


NAZZA 9G(A=2 7-1) = T = ( a) 


NA=1,Z). GA=1,F) Act kT 


where Qə is the binding energy of this new species. If you multiply 
these two equations, then you obtain a relation between A, Z, and 
A—2,Z—1 


N(A—2, Z—1)Nn Np _ZŽG(A—2, Z—1) 62 (4-2 ae 
N(A, Z) 7 G(A, Z) A 


sD (-S (26.2.4) 


This progression can be taken all the way to Z = A = 1, with the 
result 


GAZ 
QA 


N(A, Z) = > ANI NA O14 exp (-3) (26.2.5) 


kT 
where 
Q=(Zmy+(A-—Z)M, — M(A,Z))c? (26.2.6) 


is the binding energy of the nucleus A, Z. Note the resultant abun- 
dances are a function only of the temperature, the binding energy of 
the species, and the number density of protons and neutrons. Note 
also that the latter two quantities are constrained. Since the density 
of the star is known, we have from (5.1.1) 


Nap=Np+Nnt >) NA; (26.2.7) 





Moreover, if we parameterize the system with the quantity Z/N, 
t.e., the average proton to neutron ratio, then 


Z == ~-SLZN(A,Z) +N, 
N YS (A-Z)N(A,Z)+ Np 
Thus (26.2.5), (26.2.7), and (26.2.8) form a set of three equations 


with three unknowns, and p, T, and Z/N described the entire set 
of nuclear abundances. 


(26.2.8) 





Unfortunately, the analysis above is not completely accurate, since 
it fails to account for reactions that involve the weak nuclear force, 
i.e., beta decays. Since the neutrinos from these decays leave the 
star, there is no inverse reaction associated with beta-decays, and 
no equilibrium condition. To first order, however, the effect of beta- 
decays can be modeled through their impact on Z/N. Thus, 


(F) > ($) + +a? ($). (26.2.9) 


The computation of beta-decay rates is a bit tricker than normal, 
since many of the decays come from excited nuclear states, and 
the high densities increase the probability of electron captures by a 
factor of ~ 100 over terrestrial conditions. However, these reactions 
are important, in that a small change in Z/N greatly affects the 
composition of the iron group. Clearly, the faster the reactions 
proceed, the less time there is for beta-decay, and the closer Z/N 
will be to 1. 


The question of beta-decays and the value of Z/N has important 
consequences for supernova calculations and for cosmic abundance 
determinations. For example, consider silicon burning, where 7°Si 
is converted to iron peak elements. Let’s examine the case where 
the two dominant resultant species are °*Fe and °°Ni. These nuclei 
are in statistical equilibrium via the reactions 
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If ?8Si burns to Ni, then the reaction is exothermic with an net 
energy release of 10.9 MeV and the star can continue burning; on 
the other hand, if the result is °*Fe + 2p, then —1.3 MeV is lost 
in the reaction, and the star will collapse. From (26.2.4), the ratio 
between °4Fe and °°Ni is 


N(Fe) 19 _ 42 G(4Fe) a Tes a + = 


N(®6Ni) ? G(56Ni) (56 kT 


The ratios of the partition functions for these two species is ~ 1, so 
plugging in the numbers yields 


N(°*Fe) 


N2 = T3 1068-13—62.09/Tə 26.2.10 
NOON: Boe l ) 


We can now substitute for N, using (26.2.8); if Z/N =1 


Z 26.N (54Fe) + 28N (8NI) + N. 
= | a E a SN ON 
N 28N (Fe) + 28N (Ni) p ooo) 
Substituting this in (26.2.10) gives 


N(°4Fe)? _ N2 X (°4Fe)? A(°°Ni) p? = TS 1067-53-62-09/To 
N(°®Ni) X (°°Ni) A(>4Fe)3 2 
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The graph for Z/N = 1 shows that at low temperatures, 7°Si will 
fuse to °°Ni (which then beta-decays to °°Fe); this reaction pro- 
vides energy to supports the star. However, if 7°Si burns at a high 
temperature, °*Fe will result, and the star will collapse. 


This process of nucleosynthesis is call the e-process, for equilibrium 
process. 


In general, the peak iron element will have approximately the same 
proton to neutron ratio as given by Z/N. Thus, 


26 
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In nature, the most abundant isotope in the iron group is °°Fe, 
with 26 protons and 30 neutrons. This suggests that the e-process 
usually either occurs very rapidly at low temperatures, so that °°Fe 
is produced via the decay of °°Ni, or proceeds very slowly so that 
beta-decays make Z/N 7% 0.87. 
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The binding energy per nucleon of the most stable nucleus of each 
atomic weight. The solid circles represent nuclei with even numbers 
of protons and an even number of neutrons, while the crosses show 


odd-A nuclei. 
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The early phase of nuclear rearrangement of initially pure 7°Si at 
T = 5 x 10° K and p = 1.3 x 10’ g cm~’. The abundances grow 
very rapidly at first, as the liberated alpha particles are consumed 
in the rapid flow towards the iron group. 
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Number densities of key nuclei during silicon burning at T = 5 x 
10° K at later times. The densities of free protons and neutrons 
have become nearly quasistatic as a result of the near equilibrium 
between their rate of capture and their photoejection into the nu- 
clei heavier than 7°Si. All the heavier nuclei have approached a 
quasi-equilibrium with 7°Si and the pool of free nucleons, so their 
abundances also change very little. 
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Mass fractions of the most abundant nuclear species in nuclear sta- 
tistical equilibrium at T = 4 x 10° K and p = 10° g cm7?. As the 
value of Z/N is reduced, the equilibrium shifts to more neutron-rich 
nuclei (Clifford & Taylor 1965, Mem. Roy. Astron. Soc., 69, 21.) 


5 Ni56 
Fe 6 Fe54 Fe 44 2p 


log X 


2.0 3.0 4.0 5.0 6.0 
Tam 
<— Burning time 

A schematic of the principal nuclei that result from the burning of 
silicon. For T > 4 x 10° K, the burning is so fast that Z/N = 1 
throughout the burning. The time required for the burning increases 
greatly as the temperature is decreased; this results in lower values 
for Z/N. Near T ~ 3 x 10° K, the neutron-rich nucleus °*Fe has 
sufficient time to appear; at even lower temperatures the beta decays 
drive the equilibrium to °°Fe. 








